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The spin diffusion/transport in n-type (001) GaAs quantum well at high temperatures (> 120 
K) is studied by setting up and numerically solving the kinetic spin Bloch equations together with 
the Poisson equation self-consistently. All the scattering, especially the electron-electron Coulomb 
scattering, is explicitly included and solved in the theory. This enables us to study the system 
far away from the equilibrium, such as the hot-electron effect induced by the external electric 
field parallel to the quantum well. We find that the spin polarization/coherence oscillates along 
the transport direction even when there is no external magnetic field. We show that when the 
scattering is strong enough, electron spins with different momentums oscillate in the same phase 
which leads to equal transversal spin injection length and ensemble transversal injection length. It 
is also shown that the intrinsic scattering is already strong enough for such a phenomena. The 
oscillation period is almost independent on the external electric field which is in agreement with the 
latest experiment in bulk system at very low temperature [Europhys. Lett. 75, 597 (2006)]. The spin 
relaxation/dephasing along the diffusion/transport can be well understood by the inhomogeneous 
broadening, which is caused by the momentum-dependent diffusion and the spin-orbit coupling, and 
the scattering. The scattering, temperature, quantum well width and external magnetic/electric field 
dependence of the spin diffusion is studied in detail. 

PACS numbers: 72.25.Rb, 72.25.Dc, 72.20.Ht, 67.57.Lm 



I. INTRODUCTION 

The study of semiconductor spintronics 1 has caused 
a lot of attention since long spin lifetime and large 
spin transport distance have been observed in n-type 
semiconductors, 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 due to the great potential of 
the device applications such as quantum memory devices, 
spin valves and spin transistors. High spin injection 
efficiency and long spin diffusion/transport distance are 
two prerequisites of realizing these devices and have been 
widely studied both experimentally 9,10 ' 11,12 ' 13,14 ' 15,16 and 
theoretically 17 . 1 8, 19,20, 21,22,23,24,25,26,27,28,29,30,31,32,33,34,351 

Many works are focused on improving the spin injection 
efficiency by electrical method. 12,13 ' 14,30,36,37 Others 
are concentrated on the spin diffusion/transport in- 
side the semiconductors without taking care of the 
injection. 15 ' 16,20,21 ' 23 ' 24 ' 25 ' 26 ' 28 ' 32 ' 33 In the latter case, 
the spin polarization can be prepared by optical 
excitations. 15 ' 16,38 A number of theories have been 
developed to study the spin transport, such as the 
two-component drift-diffusion model, 20,22,23 ' 24 the ki- 
netic spin Bloch equation approach, 21,25 ' 27,28 ' 32 the 
Monte-Carlo simulation, 29 ' 30,31 ' 33,40 and the microscopic 
semiclassical approach. 41 ' 42,43 ' 44 In these theories, Weng 
and Wu showed that the correlations between the spin- 
up and -down states, i.e., the off-diagonal terms of the 
density matrix in spin space, play an essential role in spin 
diffusion/transport. 21 By constructing the kinetic spin 
Bloch equations by means of the non-equilibrium Green 
function method and numerically solving these equations 
with the scattering explicitly included, they predicted 



spin oscillations along the spin diffusion in the absence 

of the magnetic field 25,27 which cannot be obtained 

from the two-component drift-diffusion model. These 

oscillations were later observed in experiments. 15 ' 16 Now 

most of the theoretical works include the off-diagonal 
terms 16 , 21, 25, 26, 27,28, 29, 32, 33,41 

Another important consequence of the off-diagonal 
terms is that they allow spin to precess along the effective 
magnetic field which origins from the spin-orbit coupling 
(SOC) and is momentum dependent. The fact that spin 
with different momentum precesses with different fre- 
quency is referred to as inhomogeneous broadening. 45 ' 51 
It was shown by Wu et al. in the spacial uniform sys- 
tems that in the presence of the inhomogeneous broad- 
ening, any spin conserving scattering, including the 
Coulomb scattering, can cause irreversible spin relax- 
ation/dephasing (R/D). 45 ' 46 ' 47 ' 48 ' 49 ' 50 Later Weng and 
Wu extended this concept to the spacial inhomogeneous 
systems and showed that the spin with different momen- 
tum precesses with different frequencies during the diffu- 
sion along the spacial gradient, thanks to the off-diagonal 
term. 21 This serves as additional inhomogeneous broad- 
ening and can cause additional spin R/D combined with 
the spin conserving scattering. 21,25 ' 28 

In our previous works, we have studied both the 
transient 25 ' 28 and the steady-state 21,25,27 spin diffu- 
sion/transport in GaAs quantum well (QW) using the 
kinetic spin Bloch equations, first without 21,25 and then 
with 27 ' 28 the Coulomb scattering. The other scattering 
such as electron-phonon and electron-impurity scattering 
is explicitly included. By solving the kinetic spin Bloch 
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equations combined with the Poisson equation, one is 
able to obtain the mobility, diffusion length and spin in- 
jection length without any fitting parameter. Moreover, 
this approach is valid not only for systems near the equi- 
librium, but also for those far away from the equilibrium 
such as systems under strong external electric field (hot- 
electron effect) and/or with large spin polarization. It is 
also applicable to systems in both the strong scattering 
regime and the weak one. 50 Nevertheless, in our previous 
investigation of the steady-state spin transport, 21,25 the 
boundary conditions we used are the single-side ones, i.e., 
Pk x aa' ( x = 0) are given from the single side of the sample 
regardless of the sign of k x , with pk maa i{x) representing 
the density matrix elements. This is an approximation 
as in principal p^^ aa i{x — 0) for k x < (k x > 0) can 
only be determined from the right (left) side of x = 0. 
Moreover, by using this boundary conditions, the grid of 
the space has to be very small due to the reason spec- 
ified in the next section. Therefore the length of the 
sample we studied before is very short, even less than 
one oscillation period, as the numerical calculation is too 
time-consuming. 21,25 It is therefore also almost impos- 
sible to discuss situations such as spin transport under 
strong external electric fields. In this paper, we adopt the 
double-side boundary conditions which are widely used 



in the charge transport, 



52,53,54 



i.e., using Pk x aa'(x = 0) 



and p_fe x(T(T /(x = L) as boundary conditions with k x > 
and L standing for the right side of the sample. We 
further develop new numerical scheme to solve the ki- 
netic equations with all the scattering included. This 
allows us to solve the spin diffusion/transport with fast 
speed and high accuracy. We therefore study the spin 
diffusion/transport with large sample scale under vari- 
ous conditions such as temperature, well width, magnetic 
field and electric field. The hot-electron effect to the spin 
transport is also discussed in detail. 

This paper is organized as follows: In Sec. II we set 
up the model and construct the kinetic spin Bloch equa- 
tions with the boundary conditions by using the Kcldysh 
Green function method. In Sec. Ill we present our main 
results of the spin diffusion/ transport with different scat- 
tering, temperature, quantum well width, external mag- 
netic field and external electric field. We conclude in Sec. 
IV. In Appendix A we give the scheme for the numerical 
calculation. 



II. MODEL AND KINETIC SPIN BLOCH 
EQUATIONS 

We begin our study from a two dimensional electron 
gas (2DEG) confined by an infinite square well with 
width a along the (001) direction. The growth direction 
is assumed to be along the z-axis. A moderate magnetic 
field B and an electric field E are applied along the x-y 
plane with the x-axis being the diffusion/transport di- 
rection. The electron state is described by a subband 
index n which comes from the confinement of the QW, 



a 2D momentum k = (k x , k y ), which represents the mo- 
mentum along the x-y-plane and the spin index a, which 
stands for the spin-up (|) or -down (|) state along the 
z-axis. In the present paper, the width of the QW is 
taken to be so small that only the lowest subband n = 1 
is occupied and the higher subbands are negligible. 

By using the noncquillibrium Green function method 
with the gradient expression as well as the generalized 
Kadanoff-Baym ansatz, 55 we construct the kinetic spin 
Bloch equations as follows 21,25 



|p k (r,t) = |p k (r,i) 



• (1) 



Here p k (r, t) = ( ^e^] are the density matrices of 
VPkiT hi J 

electrons with momentum k at position r = (x,y) and 
time t. The diagonal elements h,a stand for the electron 
distribution functions of spin a whereas the off-diagonal 
elements /?ku = P^ii = Pk represent the correlations 
between the spin-up and -down states. 



dpk(r,t) 



dt 



= i{V r e k (r,i),V k p k (r,i)} (2) 



are the driving terms from the external electric field. 
Here 

k 2 <r 

e k (r,t) = _ + [ 9MB B+h(k)]---e*(r)+£ H F(r,i) (3) 
2m I 

with m* denoting the effective mass. h(k) is the 
D'yakonov and Perel' (DP) 56 effective magnetic field 
from the Dresselhaus 57 and the Rashba 58 terms. For 
GaAs, the Dresselhaus term is the leading term and h(k) 
reads 

h(k) = ( 7 Mfc 2 - <fc z 2 )),7M<fc z 2 ) - *£),o) • (4) 

Here the spin splitting parameter 7 is given by 38 



7 = {A/?>){m*/m cv ){l/^2m*3E g ){ V /^T^ri) , (5) 

in which 77 = A/(E g + A); E g denotes the band gap; A 
represents the spin-orbit splitting of the valence band; 
to* standing for the electron mass in GaAs; and m cv 
is a constant close in magnitude to the free electron 
mass to . 39 For GaAs 70 = 11.4 eV-A 3 . (fc 2 ) stands 
for the average of the operator —(d/dz) 2 over the elec- 
tronic state of the lowest subband and being (n/a) 2 
under the infinite-well-depth assumption. £hf(i", i) — 
— ^ q V r q p k _ q (r, t) is the Hartree-Fock term, i.e. the ex- 
change term, 55 with V q = £ 9z K0 (q*Zll +K *) l J (^)| 2 thc 
screened coulomb potential and k 2 — 47re 2 n / '(aksT) the 
Debye-Hucke screening constant. Here n$ is the static di- 
electric constant, n represents the 2D electron density 
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and |/(ig z )| 2 = tt 4 sin 2 y/[y 2 (y 2 — tt 2 ) 2 ] is the form fac- 
tor with y = q z a/2. The Hartree-Fock term is important 
only when there is a large spin polarization. 47 For the 
small spin polarization explored in this paper, it can be 
ignored, \l/(r) is the electric potential which satisfies the 
Poisson equation: 59 



V 2 *(r) = e[n(r) - N (r)]/(aK ) 



(6) 



with n(r) standing for the electron density at position r 
and -/Vo(r) representing the background positive charge 
density. The bracket {A, B} = AB + BA in Eq. (2) is 
the anti-commutator. 



dt 



dif 



i{V k e k (r,t),V r p k (r,£)} (7) 



are the diffusion terms which also give the additional 
inhomogeneous broadening due to the broadening of 
the spin precession frequencies for different momen- 
tum k along the spacial gradient. 21 The coherent terms 



9pk(r,t) 
dt 



and the scattering terms dpk g*' t 



can be 



found in Refs. 48 and 60. 

In the present paper it is assumed that a spin polar- 
ization is injected from the left side of the sample and 
diffuse/transport along the x-axis. Its spacial distribu- 
tion along the y-axis is uniform. The kinetic spin Bloch 
equations are then simplified into 



<9p k (x,i) d^{x,t) dp k (x,t) 



dt 



dx 



dk x 



+i 



(gp B B + h{k)) ■ -,p k (x,t) 



k x dpk{x,t) 
m* dx 
dpk(x,t) 



dt 



.(8) 



Here the spin-orbit coupling dh(k)/dk x ■ er/2 in the dif- 
fusion term is very small in comparison to k x /m* and 
is ignored. The inhomogeneous broadening can be eas- 
ily seen from these equations. For the spacial uniform 
case, dpu(x, t)/dx = and the electron spin precession is 
inhomogeneously broadened only due to the momentum- 
dependent effective magnetic field B e (k) = h.{k)/gpB + 
B. With this inhomogeneous broadening, spin conserv- 
ing scattering including the Coulomb scattering 21 ' 45,50 ' 61 
results in an irreversible spin R/D. Moreover, the scatter- 
ing can also give a counter effect to the inhomogeneous 
broadening. 50 It is noted that unless the g- factor is inho- 
mogeneously broadened, 45,62 the external magnetic field 
itself cannot cause any inhomogeneous broadening. So 
there is no spin R/D when h(k) = 0. 45 For the spacial 
inhomogeneous case, besides the spin precession around 
the effective magnetic field, its diffusion into the semi- 
conductor also accompanies a spin precession with the 
precession frequency being k-dependent, i.e., the third 
term in Eq. (8). So this causes additional inhomoge- 
neous broadening. Moreover, in the steady-state spin 
diffusion/transport, the spin precession rate along the 
diffusion direction is determined by \gpB& + h(k)\/k x 
according to Eq. (8). Therefore, differing from the spa- 
cial uniform case, the external magnetic field itself alone 



can cause spin R/D along the diffusion direction due to 
the inhomogeneous broadening of the diffusion velocity 
k x /m*. 21 

The kinetic spin Bloch equations (8) together with the 
Poisson equation (6) are highly nonlinear and have to 
be solved numerically. In our previous paper, 25 this set 
of equations are solved with single-side boundary condi- 
tions. In such a method, an extra parameter, so called 
the coefficient of viscosity, 52 is introduced to keep the 
smoothness of the solution and the Poisson equation is in- 
evitable to keep the system from divergence even though 
there are no charge imbalance and external electric field. 
This can be seen by retaining the third terms in Eq. (8) 
and replacing the scattering terms by the relaxation time 
approximation: 



k x dpk(x) 

m* dx 



Pk(x) - pi 



(9) 



with p k representing the density matrices in the equilib- 
rium. It is clearly seen from this equation that pk(x) 
is divergent when x — > oo for k x < 0. In order to cir- 
cumvent this divergence, the Poisson equation has to be 
applied and the grid for the space has to be very small. 
In the present investigation we extend the conventional 
numerical scheme of solving the Boltzmann equation for 
charge transport 52 ' 53 to solve the kinetic spin Bloch equa- 
tions by using the double-side boundary conditions. This 
conditions assume a finite sample length L and are given 

by 

{p k (x = 0, t) = F(e k;E - p(0)) , for k x > , 
\ Pk (x = L,t) =F(e KE -p(L)) , forfc x <0. 

HercF(e k;E -^(x)) = [exp(e k , E -A*(a:))/(fcBr e ) + l]" :L rep- 
resent the drifted Fermi distributions in the spin space, 

.0 Pi 

Po 
po 

left side and L is assumed to be large enough so that 
the spin polarization vanishes before it reaches the right 
boundary. e k ^ E = (k — m*^ E E) 2 /2m* is energy spectrum 
drifted by the electric field. T e and p% are the hot electron 
temperature and mobility respectively which are calcu- 
lated at the spacial uniform case following the method 
addressed in our previous paper. 48 The initial electron 
distributions are given by 



with the chemical potential matrices p(0) 
and p(L) = ( ^ u " ) if the spin injection is from the 



pu{x, t = 0) = F(e k , E - p(x)) 



(11) 



with M0) = ( M T °J and p(x * 0) = ^° ^. The 

chemical potentials Mt(1) an( l Mo are determined by the 
equations 

Y.TrfafaO)] = n(x) , (12) 

k 

^Tr[p k (0,0)a z ]-Pn(0) , (13) 
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with P denoting the spin polarization at the bound- 
ary x = 0. The double-side boundary conditions are 
in good consistence with the physical consideration that 
the electron must move along its momentum direction. 
For the Poisson equation, its boundary conditions are 
set to be \I/(0) = and = EL with E representing 

the strength of the external electric field. The numerical 
scheme for the scattering terms has been laid out in detail 
in our previous paper 48 and the new numerical scheme 
for solving the kinetic spin Bloch equations is given in 
Appendix A. By using this new scheme, we further re- 
move the extra parameter — the coefficient of viscosity. 

III. NUMERICAL RESULTS 



TABLE I: Parameters used in the numerical calculations, f2 L o 
is the LO phonon energy. 

Three quantities are used to describe the spin dif- 
fusion/transport lengths in the steady state: spin- 
dependent electron density 

N a (x)=Y,hA*), (15) 

k 

the incoherently summed spin coherence 

p{x) = Y,\p*n(*)\ > (16) 

k 



We numerically solve the kinetic spin Bloch equations 
(8) together with the Poisson equation (6) iteratively 
to achieve the self-consistent solution. We include all 
the scattering, i.e., the electron-electron Coulomb, the 
electron-phonon and the electron-impurity scattering. As 
we are interested in the diffusion/transport properties 
at high temperatures (T > 120 K), for the electron- 
phonon scattering we only need to consider the electron- 
longitudinal-optical (LO) phonon scattering. 65,66 ' 67 All 
matrix elements of the interactions and the expressions 
of the scattering terms can be found in Ref. 48. The 
width of the QW and the electron density are taken to 
be a = 7.5 nm and N e = 4 x 10 11 cm~ 2 separately unless 
otherwise specified. The initial spin polarization of the 
left boundary is P — 5%. In the numerical calculation, 
we divide the momentum k-space into N x M control re- 
gions as plotted in Fig. 8 in Ref. 48 and the corresponding 
center points are 



k„, m = ^(n+-)A£^sin— ,cos— J (14) 

with n = 0, • • • , N - 1 and m = 0, • • • , M - 1. In the 
present work, we take (N,M) = (16, 18). 68 The energy 
cut is determined by the ratio between the distribution 
at the energy cut and the distribution at the zero energy 
is less than 0.1 %, which is around lOEf at T = 200 
K with Ef = 14.3 meV being the Fermi energy. The 
sample length is usually taken to be L = 10 /im, unless 
otherwise specified. This length has to be taken long 
enough so that the spin polarization effectively vanishes 
at the right boundary. The other parameters are listed 
in Table I. 69 



K 10.8 




12.9 


Old 35.4 meV 


m* 


0.067m 


g 0.44 


E g 


1.55 eV 


a 5.331 


A 


0.341 eV 



and the coherently summed spin coherence 

p'(x) = \J2p kn (x)\ . (17) 

k 

Unlike the corresponding quantities defined in the time 
domain, 70 these three quantities are defined in the spacial 
domain. The irreversible spin dephasing length L p can be 
deduced from the decay of the envelope of p; 21 - 63 - 64 the 
ensemble spin diffusion length L* is given by the decay of 
the envelope of p'; whereas the spin diffusion length Ld 
is determined from the decay of the envelope of AN — 
iY T -A x . 



A. Transient momentum-resolved spacial evolution 
of spin polarization 

We first show the spacial dependence of /k j(T for differ- 
ent momentums k at t = 5, 10 and 300 ps in Fig. 1(a) 
with the electron-electron and the electron-LO-phonon 
scattering included. T = 200 K and B = 0. Differ- 
ing from the intuition that the spin precession should 
be controlled by the effective magnetic field h(k), it is 
seen from the figure that the electron spins with differ- 
ent momentums present almost the identical precession 
period. This is understood due to the strong scattering 
which enforces different spins to precess with the same 
frequency. In order to reveal this effect, we scale all the 
scattering terms with a scaling factor A, i.e., replacing 
the scattering terms in Eq. (8) by \ ®£h&£l | gcat . First 
we consider the case without any scattering and electric 
field, i.e., A = and E = V r *(r) = 0. It is then easy 
to obtain the analytical solution of Eq. (8) in the steady 
state with the boundary conditions (10): 



■5 



Pk{x) = < 

[pk{x = L), 



by further neglecting the HF terms. This solution is 
somehow different from our previous analytical result 25 
due to the different boundary conditions, ft describes 
the injection of the electron spin with momentum k from 
both boundaries into the 2DEG. For each electron with 
k x < 0, it comes from the right boundary (x — L) where 
there is no spin polarization. So there is no spin preces- 
sion in the 2DEG. However, for electron with k x > 0, it 
comes from the left boundary where the electron spin is 
polarized. The spin polarization precesses in the 2DEG 
due to the presence of the effective magnetic field B e (k). 
Moreover it does not decay in the absence of any scatter- 
ing, ft is noted from Eq. (18) that unlike the temporal 
spin precession where the precession rate is determined 
by <7/UB|B e (k)|, the spacial spin precession rate is deter- 
mined by 

fl k = m*g^ B B e (k)/k x (19) 

which strongly depends on the velocity k x /m*. This is 
because in Eq. (8) the diffusion term is proportional to 
the velocity. This gives strong inhomogeneous broaden- 
ing along the spin diffusion. 21 

Then we introduce the scattering into the system, first 
a small one with A = 0.05 [Fig. 1(c)] and then the normal 
one with A = 1 [Fig. 1(a)], which are compared with the 
case with A = [Fig. 1(b)]. Once the scattering term 
is turned on, the electron spin presents a very different 
evolution behavior from Eq. (18). First, electrons with 
k x < [curves for different wavevector k„ !m (Eq. (14)) 
with (n, m) — (0, 7) and (5, 7)] obtain a spin polarization 
due to the scattering which scatters electrons with k x > 
to those with k x < and hence the spin polarization is 
transferred from electrons with k x > 0. Consequently 
this kind of spin polarization is also built up from the 
left edge to the right one even though k x < 0. Second, 
the scattering tends to drive the spin precession and spin 
diffusion to reach the same precession frequency and dif- 
fusion velocity respectively for different k. This can be 
seen from Fig. 1 where A/k = /kf — /k| is plotted against 
the position x for different momentum k at t = 5, 10, and 
300 ps (already in the steady state) by comparing the re- 
sults with A = (no scattering), A = 0.05 (weak scatter- 
ing) and A = 1 (normal scattering). In the computation 
we only include the electron-electron and electron-LO- 
phonon scattering. It is easy to find that when there 
is no/ weak scattering, spin precession frequencies and 
spin diffusion velocities with different momentum show 
different values. However, when the scattering becomes 
stronger, the spacial spin precession frequencies and the 
spin diffusion velocities for different momentum k tend 
to a single frequency and velocity as shown in Fig. 1(a). 



I 

The same behavior is also observed for the spin coherence 
Pktx as shown in Fig. 1(d) in the steady state (t = 300 
ps). This is also the reason that the numerical results 
are highly accurate even for small grid numbers in the 
k-space. Moreover, it shows that when the scattering 
is strong enough, there is no interference-induced decay 
and consequently L p = L*. The similar results have been 
obtained in the time domain. 70 In fact, even for the case 
N e = 4 x 10 10 cm" 2 and T = 120 K, the intrinsic scat- 
tering is already strong enough to ensure this behavior. 
Spin echo experiments are needed to verify our findings 
here. 71 

This behavior can be understood as following: With- 
out the scattering, each spin precesses with its own fre- 
quency J7k independently with the k-dependent preces- 
sion frequency referred to as the inhomogeneous broad- 
ening. Scattering tends to suppress this inhomoge- 
neous broadening and to make it a homogeneous one. 
Our result indicates that the intrinsic scattering (the 
electron-electron and electron-phonon scattering) is al- 
ready strong enough to totally suppress the inhomo- 
geneous broadening. The spin polarizations for differ- 
ent k all precess with single frequency (fik) = = 
m *l{{kl) - (fc 2 ),0,0) when B = 0. In the case we 
calculate here, O = -1.61//xm (7 = 11.4 cV-A 3 , 
(ky) ~ ^^-ksT for the non-degenerate case) which cor- 
responds to the spacial period 3.88 /urn. 

Moreover, it is seen from the figure that besides the 
counter effect of the scattering to the inhomogeneous 
broadening, it also causes irreversible spin R/D. 45 In 
the steady state, there is no spin relaxation for each k, 
as shown in Fig. 1(b) in the absence of any scattering. 
Whereas strong spin relaxation is observed in Fig. 1(a) 
in the presence of unsealed intrinsic scattering. 

Finally the scattering also affects the total spin polar- 
ization at the left boundary in the steady state. It is 
noted that the boundary conditions Eq. (10) at x = 
are only for k x > 0. The total spin polarization Eq. (13) 
is determined by the contributions from both k x > and 
k x < 0. Without scattering, the spin polarization for 
k x < is always zero [Eq. (18) and Fig. 1(b)]. Therefore, 
the total spin polarization is only one half of 5 % at the 
left boundary. The scattering transfers spin polarization 
from states of k x > 0. Therefore, the larger the scat- 
tering, the closer the total spin polarization approaches 
the boundary condition: P = 3.0 % for A = 0.05 and 
4.5 % for A = 1.0. We believe that the spin polariza- 
tion at the left boundary can be very close to 5 % for 
sufficiently strong scattering. In such a case, the spin in- 
jection with double-side boundary conditions can be ap- 
proximated by the single-side boundary conditions used 
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FIG. 1: (Color online) A/ k [(a)-(c)] and pun [(d)] vs. x for 
different k n ,m at t = 5, 10 and 300 ps (with the system being 
in the steady state at 300 ps). T = 200 K and 5 = 0. (b) no 
scattering: A = 0; (c) small scattering: A = 0.05; (a) and (d) 
normal scattering A = 1. Only the electron-electron and the 
electron-LO phonon scattering is included in the calculation. 
It is noted that although x is plotted up to 6 |Um, the sample 
length L used in the calculation is 10 /im. 



in our previous papers. 21,25 However, for small spin po- 
larizations, the spin injection properties such as the spin 
diffusion length and the spin precession period are not 
determined by the initial spin polarization at the bound- 
ary. 



B. Effect of scattering on spin diffusion 

As shown in the previous subsection that the scattering 
plays a crucial role when the electron diffuse/transport 
into the 2DEG, we now study the effect of the scattering 
on the total spin signal. In Fig. 2 the spin-resolved elec- 
tron density N a [Eq. (15)] and the incoherently summed 
spin coherence p [Eq. (17)] in the steady state are plot- 
ted against the position x by first including only the 
electron-LO-phonon scattering (red curves), then adding 
the electron-electron scattering (green curves) and finally 
adding the electron-impurity scattering (blue curves) 
with Ni = N e at T = 120 K (a) and 300 K (b). B = 
in the calculation. 




2 4 6 

x (|xm) 

FIG. 2: (Color online) Effect of the scattering on spin diffu- 
sion in the steady state at (a) T = 120 K and (b) T = 300 K. 
Red curves: with only the electron-LO-phonon (EP) scatter- 
ing; Green curves: with both the electron-electron (EE) and 
EP scattering; Blue curves: with all the scattering, i.e., the 
EE, EP and electron-impurity (EI) scattering. The impurity 
density Ni — N e . Note the scale of the incoherently summed 
spin coherence is on the right hand side of the figure. 
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X (|im) 

FIG. 3: (Color Online) 7V CT and p vs. the position x at T = 
120, 200 and 300 K with width a = 7.5 nm and L — 10 
pm. Ni = 0. Note the scale of the incoherently summed spin 
coherence is on the right hand side of the figure. 

It is seen from the figure that similar to the transient 
spin diffusion in our previous investigation, 25 ' 27 ' 28 both 
N a and p oscillate with the position even when there 
is no external magnetic field here. Besides the oscilla- 
tion, strong spin relaxation and dephasing are observed. 
Nevertheless, they present a more complicated behavior: 
When T — 120 K, the spin relaxation and dephasing al- 
ways decrease when a new scattering is added. However, 
when T = 300 K, they become faster when the electron- 
electron scattering is added to the case with only the 
electron-LO-phonon scattering but become slower when 
electron-impurity scattering is further added. Finally, 
the spin polarizations at the left boundary are 4.75 %, 
4.20 % and 3.60 % at 120 K and 4.8 %, 4.7 % and 4.6 
% at 300 K for the cases with all the scattering, with 
both the electron-electron and the electron-LO-phonon 
scattering and with only the electron-LO-phonon scat- 
tering, respectively, i.e., the spin polarization at the left 
boundary always increases with the strengthening of the 
scattering, either the strengthening is due to the addition 
of more scattering or due to the increase of temperature. 
This is consistent with the result in the previous subsec- 
tion. 

The spin precession in the steady-state spin injection 
in the absence of the magnetic field is in consistence with 
the latest experimental observation in bulk system. 15 
This kind of spacial precession is different from the spin 
precession in the time domain, 11 where the spin preces- 
sion can also be observed in the absence of any magnetic 
field in 2DEG when T < 2 K. This is because the inho- 
mogeneous broadening in the time domain is h(k) with 
its average over a homogenous electron distribution being 
(h(k)} = 0, which is too weak to sustain the scattering 
when the temperature is higher than 2 K. However, the 
spin precession in the real space during the diffusion is 
quite different. Here the inhomogeneous broadening is 
caused by l~ik = m*h(k)/k x with (fik) = Jlo ^ 0. From 
previous subsection, we know the precessions of different 



wave vector k share the same frequency fio even with 
the normal intrinsic scattering. So the total spin polar- 
ization oscillates in the space and this kind of oscillations 
exist even in the room temperature. However, in the ex- 
periments by Beck et al., 15 as (fc 2 ) = (fc 2 ) in the bulk 
system, Slo = and the effective magnetic field is pro- 
vided by strain which is very small. This is the reason in 
bulk one has to observe the spin oscillations at very low 
temperature. 

The complicated scattering dependence of the spin 
R/D in the real space is similar to that of the spin R/D 
in the time domain, 49 ' 50 where the scattering dependence 
of the spin R/D time is different in the strong and weak 
scattering regimes. 50 In the weak scattering limit, adding 
a new scattering increases the spin R/D as it provides 
more spin R/D channel in the presence of the inhomoge- 
neous broadening. Nevertheless, in the strong scattering 
limit, the counter effect of the scattering to the inho- 
mogeneous broadening is important and adding a new 
scattering leads to a weaker spin R/D. 50 Similarly, the 
spin diffusion length, which is used to represent the spin 
relaxation in the real space, is also controlled by the in- 
homogeneous broadening and the scattering, with the 
relative strength of them dividing the system into the 
weak/strong scattering regime. At 120 K, the electron- 
LO-phonon scattering is very weak and falls into the 
weak scattering regime. Therefore the diffusion length 
decreases with the scattering. But at 300 K, the scat- 
tering becomes stronger. Our results indicate the sys- 
tem is near the transition regime between the strong and 
weak scattering and therefore the changes are marginal 
by adding new scattering. 



C. Temperature and well width dependence of spin 
diffusion 

In Fig. 3 the steady-state spin-resolved electron den- 
sity N a and the incoherently summed spin coherence p 
are plotted as functions of position x at different temper- 
atures T = 120, 200 and 300 K for QW with a = 7.5 nm. 
In the calculation only the intrinsic scattering is included. 
One finds from the figure that the spin precession period 
increases with the temperature. This is because the pe- 
riod is around (fi^ 1 ) = [m*7|(fc 2 ) - (vr/a) 2 !]" 1 . With the 
increase of temperature, (fc 2 ) increases. Consequently 
the period of the spacial spin oscillation becomes longer. 

The spin diffusion length decreases with the tempera- 
ture. This is understood that with the increase of tem- 
perature, both the inhomogeneous broadening and the 
scattering increase. In the weak scattering regime, this 
increase leads to a faster spin R/D. Nevertheless, the 
temperature dependence of the spin diffusion length is 
very mild compared to the temperature dependence of 
the spin precession period. 

We further study the well- width dependence of the spin 
diffusion in the steady state in Fig. 4, with only the intrin- 
sic scattering included at T = 120 K. One finds from the 
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FIG. 4: (Color online) N a and p vs. the position x at different 
quantum well widths a = 7.5, 10, 12.5, and 15 nm (with 
L — 10 ,20, 25 and 30 fim respectively) at T = 120 K. The 
solid curves stand for iVj and p, and the dash-dotted ones 
represent N±. Note the scale of the incoherently summed 
spin coherence is on the right hand side of the figure. 



figure that both the spin oscillation period and the spin 
diffusion length increase markedly with the well width. 
This is consistent with the spin R/D time in the time 
domain 47 and is understood due to the decrease of the 
Dresselhaus spin-orbit coupling Eq. (4) in which the lin- 
ear term is proportional to 1/a 2 . 



D. Magnetic field dependence of spin diffusion 

We investigate the effect of magnetic field on the spin 
R/D in the steady state spin diffusion. It has been 
pointed out by Weng and Wu that the magnetic field 
plays a different role in spin diffusion/transport from 
the time-domain spin precession in the spacial uniform 
case. 21 In the spin precession in the time domain, the 
magnetic field alone without the DP term will not cause 
any spin R/D due to the absence of inhomogeneous 
broadening. 45 Whereas in spin diffusion/transport, as the 
inhomogeneous broadening is caused by which reads 
m* ' g\isB /k x in the absence of the DP term and still pro- 
vides an inhomogeneous broadening leading to the spin 
R/D. 21 Furthermore, from Eq. (19) one finds that the 
inhomogeneous broadening caused by the applied mag- 
netic field depends on the direction of the magnetic field 
when it is combined with the inhomogeneous broaden- 
ing caused by the effective magnetic field m*h(k) / 'k x . 
Therefore different direction of the magnetic field yields 
different spin diffusion length. 

We reveal the magnetic field effect on the steady-state 
spin diffusion by plotting in Fig. 5 the spin diffusion 
length Ld and the spin oscillation period Lq, fitted from 
the calculated spin polarization along the z-direction by 
AN z (x)/N e = Cexp(—x/Ld)cos(2Trx/Lo + 0), as func- 
tion of the applied magnetic field at T — 120 K. The 
direction of the magnetic field is either vertical (along 
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FIG. 5: Inverse of the period of the spin oscillation L^ 1 
(curves with •) and the spin diffusion length Ld (curves with 
■ ) vs. the external magnetic field B which is applied either 
vertical (||y, dashed curves) or parallel (\\x, solid curves) to 
the diffusion direction. T = 120 K and N t = 0. The dashed 
curves are for the vertical magnetic field and the solid curves 
are for the parallel one. Note the scale of the diffusion length 
Ld is on the right hand side of the figure. 



the j/-axis) or parallel (along the x-axis) to the diffusion 
direction. 

One finds from the figure that the magnetic field leads 
to additional spin R/D and therefore the spin diffusion 
length decreases with the magnetic field, regardless of 
the direction of the magnetic field. However, it is inter- 
esting to see that the spin diffusion length has different 
symmetry when the magnetic field is vertical or paral- 
lel to the spin diffusion direction: When B is parallel 
to the y-axis, Ld(B) = Ld{—B); However, when B is 
parallel to the x-axis, Ld(B) ^ Ld{—B). This can be 
understood from the fact that the density matrices from 
the kinetic spin Bloch equations (8) have the symmetry 
Pk x ,k y ,kz( B ) = Pk x ,-k y: kA~ B ) wli en B is along the y- 
axis. This symmetry is broken if B is along the x-axis. 

Differing from the magnetic field dependence of the 
spin diffusion length, the spin oscillation period Lq de- 
creases monotonically with the vertical magnetic field, 
regardless the direction of the field. However Lq in- 
creases with the parallel magnetic field when the field 
is along the diffusion direction and less than 4 T (the 
spin precession disappears when the field is larger than 
4 T) but Lq increases with the magnetic field when it is 
anti-parallel to the diffusion direction. The spacial pe- 
riod is determined by f2o + m* g /i_eB / {k x ) with SIq being 
along the x-axis. (k x ) represents the average of k x due 
to the spin gradient and can be roughly estimated by 
(k x ) = ^2kk x Afk/J2 k Afk, which is a positive value 
due to the presence of the spin gradient. For the vertical 
magnetic field, Jlo and H/(k x ) are perpendicular to each 
other, so the mag nitude of the spin oscilla tion period is 
determined by [\J [m* ^/i^B) 2 / (fc x ) 2 + ^q] _1 wn ich al- 
ways decreases with the magnetic field. However, for 
the parallel one, these two vectors are in the same direc- 
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tion and the period is determined by \"f{{ky) — (fcf)) + 
gfiBB I (k x ) | . So the period increases with the magnetic 
field. Moreover, it is interesting to see that when B ~ 4 
T, l7((fcy} — {k1))+g^B I (k x )\ ~ and further increasing 
of the magnetic field does not lead to any spin oscillation 
due to the too large inhomogeneous broadening caused 
by the magnetic field. 



E. Electric field dependence of spin diffusion 

We now turn to study the electric field effect on spin 
transport. An external electric field is applied along the 
x-axis. When the electric field is large enough, the spin 
transport is then in the hot-electron regime. It has been 
shown in our previous investigation 48 ' 72 in the spacial 
homogeneous case that the spin precession in the time 
domain is markedly affected by the electric field. Espe- 
cially an effective magnetic field jVd^ 2 / \a 2 m*), which is 
proportional to the electric field, is induced by the elec- 
tric field combined with the DP term due to the center- 
of-mass drift velocity = [lE (and hence k x — Vd/m*) 
driven by the electric field, with /i representing the mo- 
bility. This effective magnetic field leads to a spin preces- 
sion in the absence of any magnetic field and the preces- 
sion frequency changes with the variation of the electric 
field. 48 One may therefore expect that now in spin trans- 
port, the spacial spin oscillation period is also propor- 
tional to the electric field. This is not the case as again k x 
in the spin diffusion terms, i.e., the third terms in Eq. (8), 
plays an important role in the spin diffusion/transport 
and consequently the spacial oscillation period Sl^ 1 is in- 
dependent on k x . Therefore, in spin transport the electric 
field cannot induce an additional spacial spin precession. 
It is further expected that in the steady state, (k x ) is now 
determined by both the drift velocity and the spin gra- 
dient. Therefore, if the electric field E x > (< 0), (k x ) 
is reduced (enhanced) by the electric field and the spin 
diffusion length depends on the direction of the electric 
field. 

In Fig. 6 the spin density AN in the steady state is 
plotted against the position x at different electric fields 
E = 0.5, 0, —0.5 and —1 kV/cm with only the intrin- 
sic scattering included at T = 120 K. It is seen from 
the figure that the spin oscillation period is almost un- 
changed when the electric field varies from —I kV/cm 
to 0.5 kV/cm. This is consistent with the latest experi- 
mental observation in bulk system. 15 It is further noticed 
from the figure that the spin diffusion length is markedly 
affected by the electric field. To reveal this effect, the spin 
diffusion length is plotted as function of electric field in 
the same figure. One finds when the electric field varies 
from —I kV/cm to +1 kV/cm, the spin diffusion length 
decreases monotonically. This can be easily understood 
from the fact that when an electric field is applied along 
— x-direction (+x-direction), the drift velocity caused by 
the electric field is along the x-direction (— x-direction), 
which enhances (cancels) the velocity driven by the spin 
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FIG. 6: (Color online) Spin-resolved electron density 7V CT vs. 
position x at different electric field E — 0.5, 0, —0.5 and — 1 
kV/cm and the spin diffusion length Ld against the electric 
field E at T = 120 K. N t = 0. It is noted that although x is 
plotted up to 10 /im, L = 25 /im when E = —1.0 and —0.75 
kV/cm; 20 um when E = -0.5 and -0.25 kV/cm; 10 um 
when E = 0; and 5 um when E = 0.5, 0.75 and 1.0 kV/cm. 
Note the scales of the spin diffusion length are on the top and 
the right hand side of the figure. 



gradient. Therefore, (k x ) is increased (reduced) and the 
inhomogeneous broadening is decreased (enhanced). A 
longer (shorter) spin diffusion length is then observed. 



IV. SUMMARY 

In summary, we have performed a theoretical investi- 
gation on the spin diffusion/transport in n-type GaAs 
QWs at high temperatures (> 120 K) by setting up and 
numerically solving the kinetic spin Bloch equations to- 
gether with the Poisson equation, with all the scattering 
explicitly included. A new numerical scheme is developed 
which can solve the kinetic equations with high accuracy 
and speed. 

The spin R/D mechanism in the spin diffu- 
sion/transport is due to the inhomogeneous broadening, 
combined with the spin conserving scattering. Differ- 
ing from the spin precession in the spacial uniform sys- 
tem where the inhomogeneous broadening is caused only 
by the momentum-dependent effective magnetic field due 
to the SOC, in spin diffusion/ transport, the inhomoge- 
neous broadening is caused additionally by the momen- 
tum dependence in the spin diffusion rate and therefore 
even the magnetic field alone can lead to spin R/D. Due 
to the joint effects from the momentum-dependence in 
the diffusion rate and the momentum dependence in the 
Dresselhaus effective magnetic field, there are spin oscil- 
lations along the spin diffusion even in the absence of 
the external magnetic field. The period of these spacial 
spin oscillations can be affected by the temperature, the 
well width and the electron density, but is independent 
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on the external electric field. It increases when an ad- 
ditional scattering is added or the temperature/the well 
width is increased. Moreover, it is shown that when the 
scattering is strong enough, the spin oscillations of elec- 
trons at different momentum k show the same period. 
In fact, the intrinsic scattering, i.e., the electron-electron 
and the electron-LO-phonon scattering is already strong 
enough to lead to this effect. This indicates that during 
the spin diffusion, there is no interference-induced spin 
dephasing and L p = L*. 

It is shown that the Coulomb scattering makes 
marked contribution to the spin R/D in the spin dif- 
fusion/transport. It is especially important in the hot- 
electron case as it provides the important mechanism 
of thermalization. The scattering, the temperature, the 
QW width, the magnetic field and the electric field de- 
pendence of the spin diffusion length and the spacial spin 
oscillation period are explored in detail. It is shown that 
in the weak scattering regime, the spin diffusion length 
decreases if a new scattering is added into the system, 
while in the strong scattering regime, it increases. In the 
temperature regime we are interested in, the spin diffu- 
sion length decreases with the temperature. Moreover, it 
increases with the increase of the QW width as the inho- 
mogeneous broadening decreases with the well width but 
decreases with the external magnetic field as the later 
adds a new inhomogeneous broadening into the system. 
The external electric field can effectively prolong or sup- 
press the spin diffusion length depending on whether the 
electric field is antiparallel or parallel to the spin diffusion 
direction. We believe this investigation is useful for the 



understanding of the spin transport as well as the design 
of the spintronic device. 
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APPENDIX A: NUMERICAL SCHEME FOR 
KINETIC SPIN BLOCH EQUATIONS 

The kinetic spin Bloch equations Eq. (8) can be rewrit- 
ten as: 

d Jc d 

— Pk ( x ,t) + -£— pk(x,t) =Gk[x,t] (Al) 



with G k [x,t) = + 

dr 

The discretization of the time and spatial derivatives in 
Eq. (Al) reads 



dpk(x.t) 
dt 



dpkjx.i) 
dt 



Pk(x, t + At) - p kjx, t) p k (x + Ax, t + At) - p k (x + Ax, t) 



At 

^{G k [x,t]+G k [x + Ax,t\) . 



At 



k x p k (x + Ax,t + At) - pk(x, t + At) 



TO* 



Ax 



(A2) 



For k x > 0, the electron is propagated from its left positions to the right ones. Therefore, from Eq. (A2), the iterative 
format is 

1 - r 1 At 

Pk {x + Ax,t + At) = --—p k {x,t + At) + -—(p k (x,t)+p k (x + Ax,t)) + -—(G k [x, t] + G k [x + Ax, t}) .(A3) 
1 + r 1 + r 1 + r 

For k x < 0, the electron state is propagated from its right positions and the iterative format is then 

1 — r 1 At 

p k (x,t + At) = -— — p k (x + Ax,t + At) + -—(p k (x + Ax,t)+p k (x,t) + — — (G k [x, t] + G k [x + Ax, t]) .(A4) 
1+r 1+r 1+r 



r 



In these equations r = 2 l fc »l At . 

m* Ax 



the local source term Gk[x, t}. With the combined effects 
of the spin-flip terms and the scattering terms, the spin 
signal decays when it diffuses into the QW in the scale 
of the diffusion length L4. When L ^> Lj, one can study 



We truncate the spacial variable x from to L. At 
the left edge, the boundary conditions are given only for 

the states with k x > and at the right edge only for the spin injection properties, 
those with k x < 0. Note the state with k x = does 
not contribute to the diffusion process but contributes to 
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The computation is carried out in a parallel manner in points in the real space, it takes about 7 hours for the 
the "Beowulf" cluster. For a typical calculation with the system to evolute to 300 ps with a time step of 0.023 ps 
partitions 16 x 18 grid points in the (k, 6»)-spacc and 400 by 6-node AMD Athlon XP3000+CPU's. 
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